Intermediate dynamics between Newton and Langevin 
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A dynamics between Newton and Langevin formalisms is elucidated within the framework of the 
generalized Langevin equation. For thermal noise yielding a vanishing zero-frequency friction the 
corresponding non-Markovian Brownian dynamics exhibits anomalous behavior which is character- 
ized by ballistic diffusion and accelerated transport. We also investigate the role of a possible initial 
correlation between the system degrees of freedom and the heat-bath degrees of freedom for the 
asymptotic long-time behavior of the system dynamics. As two test beds we investigate (i) the 
anomalous energy relaxation of free non-Markovian Brownian motion that is driven by a harmonic 
velocity noise and (ii) the phenomenon of a net directed acceleration in noise-induced transport of 
an inertial rocking Brownian motor. 
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I. INTRODUCTION 

The phenomenon of Brownian motion has assumed a 
fundamental and influential role in the development of 
thermodynamical and statistical theories and continues 
to do so as an inspiring source for active research in var- 
ious fields of natural sciences The Brownian motion 
dynamics can conveniently be described by a generalized 
Langevin equation (GLE). The GLE was originally de- 
rived by Mori Q, Kawasaki Q, and Zwanzig [|| by use 
of the Gram-Schmidt procedure. It was further investi- 
gated by Lee using the recurrence relations method ||. 
Starting out from the well-known system-plus-oscillator- 
reservoir model as e.g. detailed in Ref. 0, @ , one obtains 
the GLE derived from first principles. The validity of a 
thermal GLE is typically restricted to the case with a 
thermal equilibrium; e.g. see in Refs. JEHM. Specif " 
ically, such a GLE dynamics reads [1, H Q : 

mv(t)+m [ j(t - t')v(t')dt' +d x U(x,t) = s(t). (1) 



Notably, the thermal noise e(i) is not correlated with the 
initial velocity, i.e., (v(0)e(t)) = 0, see in Refs. d 0] 
and in section III below. In contrast, the initial posi- 
tion x(0) typically is correlated with e(t). The memory 
friction j(t — t') is in thermal equilibrium related to the 
correlation of stationary random forces @, § ■ Kubo Q 
has addressed the common behavior of a classical equilib- 
rium bath by setting (e(i)e(t')) = mkBT^(t — t'). Here 
ks is the Boltzmann constant and T denotes the bath 
temperature. The one-sided Fourier transform of 7(i) 



'Electronic address: jdbao@bnu.edu.cn 



obeys Re7(w) > for real- valued u. This correlation re- 
sult for the thermal noise e(t) is commonly termed the 
fluctuation-dissipation theorem (FDT) of the second kind 
[8|. The nonlinear GLE can also be extended to account 
for a nonlinear system - linear bath interaction [1 H E3] , 
yielding a structure as in Eq. (1), but now with a nonlin- 
ear, coordinate-dependent friction function. It even can 
be generalized to arbitrary nonlinear system - nonlinear 
bath interactions containing then the potential of mean 
force 

With this work we aim at extending the theory of clas- 
sical Brownian motion by focussing on the intricacies 
of a possible non-Markovian with an incomplete, non- 
Stokesian dissipative dynamics [H, El El El Ell El • 
We will demonstrate that the commonly stated condi- 
tions for the equilibrium bath are generally not complete 
within the framework of linear response theory. This is 
so, because the existence of anomalous diffusion has not 
been considered in the original treatment by Kubo and 
others. Moreover, we discuss also the influence of initial 
correlation preparation between the system and the heat 
bath upon the asymptotical behavior of the force-free 
system. 



II. BIASING GENERALIZED BROWNIAN 
MOTION 

Let us first consider a free Brownian dynamics with 
U(x, t) = 0, possessing via the FDT of the second kind 
a finite-valued zero-frequency friction. If this dynam- 
ics is next subjected to a constant external force, i.e. 
—dll(x,t)/dx = F, the acceleration vanishes in the case 
of a Stokesian friction, because the external force bal- 
ances the friction force. A problem of broad interest is 
whether there exists an intermediate situation between 
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the Newtonian mechanics and such an ordinary Langevin 
formalism. This in turn necessitates a non-Stokesian dis- 
sipation mechanism such that the asymptotic long-time 
statistical probability will typically approach a station- 
ary state that explicitly depends on the initial prepara- 
tion. It is thus of great practical interest to research what 
kind of heat bath can take on this role. Such non-ergodic 
non-equilibrium thermodynamics presents a timely sub- 
ject that is presently hotly debated, both within theory 
Hi El Q El 03 HI HI and experiment H H H3 . 

For a GLE subjected to a constant force; i.e. U(x, t) — 
—Fx, the solution of (1) can be written as 



x(t) = x(0) + v(0)H(t) 



1 



1 



dt'H(t-t')(e(t') + F), 



v(t) = v(0)h(t) + — / dt'h(t-t')(e{t')+F), 



(2) 



where H(0) = and h(0) — 1. The two response func- 
tions H(t) and h(t) are the inverse of the Laplace trans- 



forms H(s) 



+ s^f(s)] and h(s) = [s + j(s)]~ 



respectively, where j(s) is the Laplace transform of mem- 
ory friction kernel, i.e. j(s) — J °° j(t) exp(—st)dt. 
Under the assumption that the characteristic equation 
s + j(s) = possesses a zero root, i.e. s = 0, the residue 
theorem implies that h(t) = b + J2i Tes[h(si)] exp(sit) 
and H(t) = c + bt + s~ 1 res[h(si)} exp(sjt). Here, Si 
(ReSi < 0) denote the non-zero roots of the above char- 
acteristic equation and res[- • •] are the residues. Within 
this context, the two relevant, generally non- vanishing 
quantities b and c are determined to read: 



1 



1 



l+7'(0)' 



^f'(O). 



(3) 



This result then requires that 7(0) = J °° j(t)dt = 0, i.e., 
implying a vanishing effective friction at zero frequency. 

The average velocity and the average displacement un- 
der the external bias F emerge at long times as 



{<«(*- 00))} = &(mo)} + ^*J + 

oo)» = {x(0)} + b({v(0)}t 



F 



1 F 
i 

2 m 



F \ F 
< \ {v(0)} + — t) + —d, 
ml m 



(4) 



(5) 



where d = — J2i res [h( s i)]/ s i is a noise-dependent quan- 
tity. Herein, we indicate by {• • •} the average with respect 
to the initial preparation of the state variables, i.e. an 
average over their initial values and (• ■ •} is the noise av- 
erage. The dissipative acceleration of a Brownian particle 
of mass m subjected to a constant force F then reads: 



F 



(6) 



where generally < b < 1 . This quantity b will be termed 
the dissipation reducing factor, the dissipation is reduced 



as b is increased. This result is intermediate between 
a purely Newtonian mechanics (obeying b — 1) and an 
ordinary Langevin dynamics (with 6 = 0) including of 
GLEs. 

The two limiting results for the asymptotic dynamics 
are found to read: (i) The Newton case with 6=1 imply- 
ing no dissipation, i.e. j(t) = with c = d — 0, yielding: 
a = F/m, v(t) = v(0) + ^4, and x(t) = x(0) + v(Q)t + 
\-^t 2 . (ii) The commonly known, ordinary Langevin sit- 
uation is obtained with 6 = 0; i.e. Eq. (3) then looses 
validity because of 7(0) ^ 0. For this case c = 7~ 1 (0) 
where 7(0) denotes the Markovian friction strength, re- 
sulting in a = 0, (v(t — > 00)) = F/(m A f(0)), and 
{(x(t - 00))} = {x(0)} + M0)}/ 7 (0) + £t/ 7 (0) + £d. 

In the unbiased case, we derive the two-time velocity 
correlation function (VCF) of free generalized Brownian 
motion in a generic form, i.e., 



{(v(ti)v(t 2 ))} = ^hQh-hD+UAo)}-^- 

h(h)h{t 2 ). 



(7) 



Here we used only the condition: (v(O)e(t)) = 0. Note 
that depending on the specific choice for the initial prepa- 
ration this velocity correlation generally is not time- 
homogeneous. The stationary velocity correlation func- 
tion becomes again only a function of \t\ — t 2 \ for the 
case that we use the equilibrium preparation with an 
initial velocity variance in accordance with the ther- 
mal equilibrium value, i.e. {w 2 (0)} = feT/m pll ]. 
The deduced asymptotic stationary VCF then reads: 
C vv (00) = 6-^p ^ 0. This causes a breakdown of the er- 
godic equilibrium state because of the initial preparation- 
dependence, which is encoded in the u(0)-dependent 
asymptotic results: (v) st = bv(0) ^ and {(v 2 )} st = 
k B T/m + b 2 {{v 2 (0)} - k B T/m). 

Likewise, the mean square displacement (MSD) of the 
force-free particle is written as 

{{x 2 (t))} = {x 2 (0)} + {v 2 (0)}H 2 (t) + 2{x(0H0)}H(t) 



m jo 
k B T 
m 



dt'H{t-t'){(x(0)e(t'))} 



2 [ H(t')dt' - H 2 {t) 
Jo 



(8) 



where the fourth term denotes the effect of initial cou- 
pling between system and heat bath. We will discuss the 
correlation preparation in the following section. Note 
that here the largest power in the temporal variation of 
the MSD involves the square of time. The averaged dis- 
placement can be related to the MSD via the generalized 
Einstein relation, reading 



K2/H2(F -> 0) = k B T eS 



(9) 



The ballistic diffusion coefficient is k 2 = 
lim t ^ 00 {(a; 2 (t))}i? = o/(2t 2 ), being related to the 
increasing rate of linear mobility ^ 2 {F — ► 0) = 



3 



\im t -,oo{(x(t))}/(Ft 2 )F=o- Here, this effective temper- 



ature formally reads: T eff := T + b(m{v 2 (0)}/k B - T), 
where T is the temperature for the common case with 
6 = 0. 

To assure the equilibrium behavior of this generalized 
Brownian motion the usual condition of Kubo's FDT of 
the second kind for the thermal noise must be comple- 
mented as follows: Consider the Fourier transform 7(0;) 
of the memory damping kernel, i.e., 



7(cj) = 7(s = -iuj), 



The real part of the former quantity is the spectral den- 
sity of noise. A ballistic diffusion with 6^0 thus requires 
that the lowest power of 7(s) is of first-order in s, imply- 
ing that the lowest power of Hej(oj) is proportional to 
lu 2 at low frequencies. Therefore, for a genuine thermal 
noise driven, force-free particle approaching at the equi- 
librium state the usual conditions must be completed by: 
lim s ^ (7(s)/s) — > 00, or lim^ (Re7(cj)/w 2 ) -» 00. 



III. INITIAL CORRELATION BETWEEN 
SYSTEM AND BATH 

Starting from the system-plus-reservoir model, one 
knows that the coupling between the system and envi- 
ronmental degrees of freedom there exist four kinds of 
coupling forms which do not involve a renormalization 
of potential or mass of the system. In these cases the 
heat bath consists of a set of independent harmonic os- 
cillators with masses rrij and oscillation frequencies u>j. 
Pervious work [2lj has shown that the random force is 
independent of the system variables for the coordinate- 
velocity coupling. Nevertheless, however, the expression 
of thermal noise will depend on the initial preparation of 
system for the coupling between the system coordinate 
(velocity) and the environmental coordinates (velocities) 
considered here. To the best of our knowledge, only a 
small number of prior studies [25L HB| have considered the 
general consequences of the detailed initial preparation 
procedure in view of the asymptotic statistical results of 
the system. 



A. The coordinate-coordinate coupling 

For a bilinear coupling between the system coordinate 
x and the heat bath's coordinates qj , the total Hamilto- 
nian can be written as: 



H=^-+U{x,t)+Y J 
2m t-r 1 



2m j 



2 3 3 



1i 



(11) 

Here and below the momenta of the system and bath's 
oscillators are related to P x = mv and pj = mjcjj, re- 
spectively, and the set Cj denote the coupling constants. 
The equation of motion of the system obeys the form of 



the GLE (l)and the thermal noise appearing in Eq. (1) 
emerges as 0, [|| [2l[ 



e(i) = -m7(t)s(0)+£b a th(*)> 



(121 



■ cos ujjt and £bath(*) * s deter- 



where 7 ft) = — - - 

J J 

mined by the initial coordinates and the velocities of the 
oscillators of the heat bath. The bath part £batri(i) °^ 
the noise explicitly reads: 



(10) £bathW=E c ^ U(0)cosc^ + 



smojjt ) . (13) 



Physically, this thermal noise e(t) obeys statistical 
properties that derive from the canonical, thermal equi- 
librium distribution of the total, combined system-plus- 
bath [1, [H, H3] : this thermal noise then again yields a van- 
ishing mean and its correlation obeys the thermal FDT 
Q. The statistical quantities involving noise the system 
variables are strictly determined by the joint probability 
(25| . The correlation involving initial position of the sys- 
tem and the thermal noise e(t'), i.e., the fourth term in 
Eq. (8) reads 



dt'H{t-t'){(x(0)e{t'))} 



= -2{x 2 (0)} / dt'H(t-t')^(t') 
Jo 

= -2H(t) * 7 W{z 2 (0)} = -2(1 - h(t)){x 2 (0)}, 

(14) 

wherein " * " denotes the convolution integral. It reads, 
* 7(*) = 35 ? /ds#(*)7(*)exp(rf) = 1 - H(t) = 
1 — h(t). Here we have used the relation below Eq. (2), 
i.e., H(s)j(s) = —sH(s). This contribution assumes 
a finite value, i.e., — 2(1 — b){x 2 (0)} in the long-time limit. 

Under the usual assumption that the thermal noise e(t) 
and the initial velocity v(0) of the system are not corre- 
lated, we find from Eq. (12) that the initial coordinate 
of the system must be uncorrelated with its initial veloc- 
ity for the coordinate-coordinate coupling case, namely, 
{x(0)v(0)} = 0, being the case for a canonical thermal 
equilibrium, cf. the Hamiltonian in Eq. (11). Therefore, 
the third term in Eq. (8) also vanishes. In the follow- 
ing we shall not consider preparations with such initial 
correlations between the initial coordinate x(0) and the 
initial velocity v(0). 



B. The velocity- velocity coupling 

For a bi-linear coupling between the system velocity 
and the velocities of the bath oscillators the total Hamil- 
tonian reads 



H 



P 2 
2m 



-u(x,t)+£ 



1 

2m a 



Pj 



m 



+ ^mj^q 2 



(15) 
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where dj denote the corresponding the coupling con- 
stants. We can derive again the GLE (1) describing the 
motion of the system with the thermal noise term now 
given by 



s(t) = v(0)J2 

3 

+ 



sin ujjt 



,t+ -^sinw,f 



qj (0) cos ojjt + — ' ' sin c 



u(0)m / l(t')dt' + Z hath (t), 



(16) 



This result evidences that the system can not arrive at 
the equilibrium state for any initial preparation of the 
particle velocity if b ^ 0. Therefore, for the validity 
of FDT one has to use at initial time a preparation of 
thermal equilibrium for the system and the heat bath. 
Nevertheless, one needs not to worry that the noise is 
uncorrelated with the initial velocity of the system for 
a common non-Markovian dynamics with b = 0. Then, 
the FDT is valid independent of the coupling form be- 
tween system and bath whenever the effective Markovian 
damping of the system is finite at zero frequency. This 
is so because the first and the third term in Eq. (19) 
vanishes for 6 = 0. 



where 7ft) = — V ■ - J ^ J - cos cod and in addition we have: 

{<*(°)£bath(<)>} = {M°)£bathW>} = °- We rec l uire 
that the FDT of the second kind is obeyed, namely 

that (e(t)) = and (e(t)e(t')} = k B Tj(t - t'). This 
is guaranteed when (qj(0)) = (4j(0)) = (^(0)^(0)) = 0, 
(©(0)^(0)) = b^Sij, and (I (0)^.(0)) = where 

1(0) = ft(0) - d//mjv(0) & 

In the case of velocity-velocity coupling, the fourth 
term in Eq. (8) vanishes if again {x(0)v(Q)} — 0. An 
additional term emerges, however, for the mean squared 
displacement (MSD) of the force-free particle due to the 
thermal noise e(t) which now depends on the initial par- 
ticle velocity v(0). From Eq. (2), we obtain 



{<* 2 (*)»add = -#(*)/ dt'H(t-t'){(v(0)s(t'))} 

Jo 

= 2H(t){v 2 (0)} [ dt'H{t-t') [ dwy{u) 
Jo Jo 

= 2(t-H(t))H(t){v 2 (0)}. (17) 

Indeed, the ballistic diffusion arises also in this case. No- 
tably, an additional contribution to the mean square ve- 
locity of the force-free particle [c.f. Eq. (7) at t\ = £2] 
emerges in this case. It reads 

{(« 2 W)}add = ZlW) f 'dt'h(t-t'){(v(0)e(t'))} 



2{v 2 (0)}h(t) / dt'h(t-t') / j(u)du 



(18) 



2h(t)(l-h(t)){v 2 (0)}, 



where we have used the inverse Laplace transform of the 
convolution integral in Eq. (18) by making use of the 
relation, h(s) A f(s) = 1 — sh(s). 

In particular, the mean squared velocity of the force- 
free particle emerges in the long-time limit as 



{{t^t ->«>)>} = 26(l-6){« 3 (0)} + 



{w 2 (0)} _MT)/r 



(19) 



IV. TEST BED FOR NON-STOKESIAN 
DISSIPATIVE DYNAMICS 

Given the FDT of the second kind by Kubo we investi- 
gate next unbiased, non-Markovian Brownian motion in 
(1) that is driven by colored noise known as the harmonic 
velocity noise (HVN) e(t) [28|, which, however, does not 
obey the above additional requirements. The HVN itself 
is produced from a linear Langevin equation, namely, 



y = e, e = -Te-n 2 y + ^t), 



(20) 



where £(t) denotes Gaussian white noise of vanishing 
mean with (£(t)£(t')> = 2riT 2 k B TS{t - t'). The coef- 
ficient rj denotes the damping coefficient of the system 
corresponding to the thermal white noise; T and £1 de- 
note the damping and the frequency parameters. The 
second moments and the cross- variance of y(0) and e(0) 
obey: {y 2 (0)} = T^fT^T, {e 2 (0)} = r)Tk B T, and 
{y(0)e(0)} = 0. The Laplace transformation of the mem- 
ory damping kernel reads 7(s) = r]Ts/(s 2 +Ts + fl 2 ) with 



Re7(w) 



r]T 2 iv 2 



{n 2 -lu 2 ) 2 + r\ 



(21) 



respectively. The latter corresponds to the spectrum of 
HVN which indeed vanishes identically at zero-frequency. 
In this case the dissipation reducing factor emerges as 
b = (1 + ^m^y 1 , and likewise, c = (1 - b) 2 /rj, d = 

(i-6)^r)- l -(i-6) ? - 2 ]. 

Using models with a bi-linear coordinate system-bath 
coupling the dynamics can be characterized by the spec- 
tral density of bath modes, J (a;), being related to 

Rej(co) = J(w)/mw [H, III Hi Si- Thus, for a weak 
coupling to a bath, as it can be realized either with 
optical-like bath modes [12, HH , broadband colored noise 
|33j | , or also for the celebrated case of a black-body radi- 
ation field of the Rayleigh- Jeans type 34] the static fric- 
tion vanishes. Yet other physical situations that come 
to mind involve the vortex diffusion in magnetic fields 
[HI , or open system dynamics with a velocity-dependent 
system-bath coupling [II HI HE HI- 

The non-Markovian Brownian motion can equivalently 
be recast as an embedded, higher-dimensional Markovian 
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FIG. 1: (color online). Behavior of diverse mean energies (see 
text) for a nonergodic, force-free Brownian particle vs. time 
t. The parameters used are m — 1.0, ksT — 1.0, r\ = 0.2, 
T = 5.0, n = 1.0, x{0) = 0, and (i) v(0) = y/2 and (ii) 1.0, 
respectively, for the total energies given by the solid lines and 
the remnant energies by the dashed lines, from top to bottom. 
The open circles is the adsorbed energy from the bath for all 
cases and also the total energy for the v(0) = case. 



FIG. 2: (color online). Accelerated, (time and ensemble)- 
averaged Brownian motor velocity in a rocking ratchet that 
is driven by HVN. The used parameters are x(0) = v(0) = 0, 
m = 1.0, k B T = 0.5, r) = 5.0, T = 22.0, Q 2 = 40.0, t p = 
25.0, A = 2.0, 4.0, 6.0, 10.0, and 15.0, from top to bottom. 
The particles undergo a finite acceleration a (times 100) vs. 
temperature T, being depicted with the inset for Aq = 10.0. 



process. The Fokker-Planck equation (FPE) for the prob- 
ability density P{x, v, w, u, y, s; t) corresponds then to the 
dynamics of a set of coupled Markovian LEs involving the 
auxiliary-variables (w,u,y, e): It obeys dtP — LppP, 
where Lpp is the associated FPE operator; when being 
formally supplemented here with a nonvanishing poten- 
tial U(x,t) 7 it explicitly reads: 

Q Q 

Lfp = -v— -m- 1 (-U'(x,t)+w) — 
ox ov 

d 

+(Tw + rjTv + fl 2 y + u)—— 



dw 



d 



d 



^ {w - £) d-u +£ dy- 



d d 2 d 2 

+ (Te + n 2 y)- + vT 2 k B T— + ^ 2 k B T—. 

(22) 



A. The energy relaxation 



- m {(v(t)) 2 } + - m {((v(t)-(v(t))) 2 )}. 



(23) 



The first part describes the remnant initial kinetic en- 
ergy of the particle, being dissipated partly by the heat 
bath environment. This part vanishes in the ordinary 
case with b = 0. The second part denotes the energy 
provided from the heat bath. It is independent of the 
initial particle velocity, but does not relax, however, to- 
wards equilibrium. The absorbed power of the particle 
from the heat bath, namely, the rate of work being done 
by the fluctuation force [38] , is 



^abs 



= lim 



7 (t-t'){{v(t)v(t'))}dt' 
i]k B T 



i + ? / /(2r) + 20 2 /( ?? r) ' 



(24) 



We use this form for the analysis of the energy relax- 
ation. The mean total energy of the thermal HVN-driven 
force-free particle reads 

{(E(t)}} = \m{{v 2 {t))} 



Note that it falls short of the equilibrium value -qkBT. All 
quantities plotted are dimensionless. Our numerical re- 
sults are depicted with Fig. 1. These results are obtained 
via the simulation of a set of Markovian LE's which are 
equivalent to the FPE in (22) with U(x, t) = 0. 
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B. Brownian motor exhibiting accelerated 
transport 

A most intriguing situation refers to Brownian motors 
[3^ | when driven by Brownian motion that exhibits bal- 
listic diffusion. Take the case of thermal HVN driving a 
Brownian particle according to (1) with a periodic po- 
tential that breaks reflection symmetry, namely 

U(x,t) — C/ [sin(27ra;)+ci sin(47rx)-|-C2 sm(6nx)]+A(t)x, 

(25) 

where U = 0.461, a = 0.245, and c 2 = 0.04 gpj. A(t) is 
a square-wave periodic driving force that switches forth 
and back among A(t) = Aq when 2nt p < t < (2n + l)t p 
and A(t) = -A when (2n + l)t p < t < 2(n + l)t p . The 
key challenge is whether a non- vanishing non-equilibrium 
current emerges that can be put to a constructive use in 
order to direct, separate or shuttle particles efficiently 
& 

Figure 2 depicts the ensemble- and driving-phase- 
averaged velocity (with the latter being equivalent to 
an average over the temporal driving period, see in Ref. 
EH), i.e., 

/t-\-2t p 
dt'{(v(t'))} (26) 

for various strengths of the driving force A , as obtained 
via the simulation of the Markovian LE's correspond- 
ing to the equivalent higher-dimensional FPE in (22). 
A startling finding is now that this averaged velocity is 
no longer a constant but rather increases linearly with 
time. This is in clear contrast to the behavior of or- 
dinary Brownian motors [39j . This directed accelera- 
tion is presented by the slope (see dashed lines in Fig. 
2) of the average velocity. For a weak rocking (i. e. 
small Aq) the phenomenon of directed motion involves 
the surmounting of barriers, thus hindering transport. In 
contrast, for a strong super-threshold rocking the aver- 
aged displacement is related to the mean square displace- 
ment via the modified Einstein relation in (9), yielding 
{{x(t))} oc — F b q t 2 in the long-time limit. Here F g g- is 
an effective tilting force stemming from the rocking of the 
ratchet potential. Amazingly, this Brownian motor can 
be accelerated because the driving and the noise-induced 
effective tilting force supersedes the acting friction force. 



V. CONCLUSIONS 

We have researched within the GLE-formalism in 
(1) an intermediate dynamics proceeding between New- 
ton and Langevin. The emerging non-equilibrium fea- 
tures are manifested by the initial preparation-dependent 
asymptotic stationary state, which is directly related 
to a non-Stokesian dissipative phenomenon which stems 
from a vanishing effective Markovian friction at zero- 
frequency. It has been found that the fluctuation- 
dissipation theory is valid and independent of the cou- 
pling form between system and bath when the effective 
Markovian damping is finite at zero frequency. In order 
to assure the equilibrium behavior of a system, the usual 
condition for the heat bath, i.e., the Kubo's fluctuation- 
dissipation theorem of the second kind, must be com- 
pleted by an additional requirement: lim s ^o(7(s)/s) - * 
00 or lim (i) _»o(R- e 7( w )/ 1 ^ 2 ) — * 00, where j(s) and 
are the Laplace and Fourier transforms of the memory 
damping. However, the system can not arrive at the 
equilibrium state for any initial preparation of the sys- 
tem if the the condition is not obeyed for the bilinear 
coupling between the system velocity and the velocities 
of environmental oscillators. 

Our findings exhibit anomalous super-diffusion in the 
form of a ballistic diffusion. Yet another riveting result is 
that the corresponding Brownian dynamics for a rocking 
Brownian motor exhibits a distinct, accelerated velocity, 
rather than the constant drift which typifies the situation 
with a Stokesian finite zero-frequency dissipation. We are 
also confident that our present results will serviceably 
impact other quantities of thermodynamic and quantum 
origin. Thus, this field is open for future studies that in 
turn may reveal further surprising findings. 
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